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ABSTRACT 


The goal of this thesis is to investigate the initial search planning phase of the Search 
and Localization Tactical Decision Aid (SALT) developed by »METRON, Incorporated 
of McLean, Virginia. SALT is a Computer Assisted Search (CAS) program intended for 
use by P3 UPDATE IV crews to assist them in optimal deployment of a sonobuoy field 
to prosecute a submarine threat. 

The initial search planning phase of SALT takes as user inputs environmental data, 
an initial elliptical Search Probability Area (SPA), an assumed target motion model, and 
the duration of the search. Outputs include a recommended sonobuoy pattern and the 
probability of detection of this pattern. The investigation of this phase of the algorithm 
is conducted in two parts, first, a scries of simulation routines is used to ensure that the 
probability of detection of the sonobuoy patterns generated by SALT is mathematically 
correct. Second, these same simulation routines arc used to determine if there arc al¬ 
ternate sonobuoy patterns that result in higher probabilities of detection. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While ever}' effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er¬ 
rors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. INTRODUCTION 


Tactical Decision Aids (TDA’s) have been developed and implemented on comput¬ 
ers in various forms to assist users in making decisions on a wide variety of topics. One 
prominent area where TDA's have proven to be of value is in the area of Antisubmarine 
Warfare (ASW). In the search theory phase of the ASW problem, in particular, there 
is an extremely large amount of data which must be assimilated and presented to the 
user in a manner conducive to rapid evaluation so that a sound decision can be made 
and acted upon. Several TDA's have been developed to aid the user in the area of search 
thcor\’. These TDA's have been grouped together under the title Computer Assisted 
Search (C.^S). A common thread which runs through all of the CAS programs is the 
concept of a "probability map". It is this probability map which allows all of the input 
data, mathematical calculations, and time step updates to be output as a 2-dimcnsional 
graphic di.eplay that is easily understood by the user in terms of relative probabilities on 
the location of the target of interest. 

Many CAS programs represent the stochastic aspect of target motion by a Monte 
Carlo simulation of a large number of possible target tracks. As information, either 
positive (detection) or negative (no detection after a period of lime), is obtained from tlie 
scaich, the probabilities of the individual tracks are updated and displaced as a proba¬ 
bility map at the user requested time. The most prominent TD.A's of this t\pe ai'c 
VPC/\S and P.ACSIi/'vRCil. This type of TD.A possesses an inherent simplicity which 
promotes ease of understanding of the underlying concepts. The use of the Theorem of 
'Iotal Probability and Bayes' Theorem arc easily applied to small scale problems which 
can be worked by hand and the results natur'ally gcnerali/e to large scale Computer As¬ 
sisted Search problems. The dinicully in this type of approach is that the specific ti'acks 
selected for the probability tuap do not cover all of the possibilities of target motion. 
[Ref. l;pp. 11-1.11-2J 

In an attempt to negate the above difficulty with the Monte Carlo approach to 
target motion, a now algorithm entitled Search and Localization Tactical Decision Aid 
(SALT) has been developed by ML! RON, Incorporated of McLean. Virginia. SALT 
is a CAS algorithm intended for use by P3 L'PD.VTL IV crews to assist them in the op¬ 
timal deployment of a soiiobuoy field to prosecute a submarine threat. The SALT al¬ 
gorithm uses a discrete time and space Markov Chain to simulate target motion. The 


1 






'tatc space of this Markov Chain is a five element vector (position, velocity, depth). The 
number of discrete cells in the state space is 21 entries each for latitude and longitude, 
eight entries each for north-south and east-west velocity components, and two entries 
for depth. Thus the current state of the target is represented probabilistically as 56,448 
five-dimensional cells. This anahtic approach to solving the problem of unknown target 
motion allows for many more possibilities of target motion than the simplistic track 
technique used in previous CAS models. [Ref. 1: pp. 11-46,11-47] 

The initial search planning phase of S.-\LT takes as user inputs environmental dam, 
an elliptical Search Probability Area (SPA), an assumed target motion model, and du¬ 
ration of search. The algorithm then outputs a recommended sonobuoy deployment 
pattern along with the probability of detection of this pattern. [Ref 2: p. 1] This initial 
search planning phase is the topic which will be investigated here. The investigation 
consists of two parts, first, the sonobuoy patterns recormnended by SALT will be tested 
using FORTR-AN simulations to ensure that the probability of detection for each pat¬ 
tern as output by S.ALT is mathematically accurate. Second, alternative sonobuoy pat¬ 
terns will be explored with the goal of developing patterns that improve on the S.ALT 
probability of detection. 

We shall begin in Chapter II by discussing the three types of motion models used 
by S.ALT to develop initial probabilits maps. Chapter Hi follows with a prescntaiion 
of how buoy patterns arc developed and esaluatcd based on the input data for each type 
of target motion. Chapter IV provides the documentation of the rORTR.'\N simu¬ 
lations that are used to verify the probabilitv of detection for the SALT generated pat¬ 
terns. In Chapter V. the results obtained from the FOR'I R.AX simulation will he 
presented and compared to the output from the S.AL’f algorithm. Chapter VI explores 
the use of alternative sonobuoy patterns with the goal of developing patterns which will 
improve the probability of detection of those patterns currently used in the SALT pro¬ 
gram. Chapter VI1 is a summary of the results. The inputs to the SALT algorithm used 
to obtain the sonobuoy patterns for testing arc listed in Appendix A. 'I he source code 
for examples of the FORTRAN simulation routines is presented in Appendix B. 





II. SALT MOTION MODELS 

There are three different models of target motion built in to the SALT algorithm. 
These models are a patrolling target, a transiting target, and a fleeing datum target. The 

specifics of each type of target motion are discussed below. In general, each scenario 

\ 

begins with an elliptical SPA which is input by the user. To specify a SPA requires seven 
input parameters: date-time-group, latitude and longitude of the SPA center, semi-major 
axis (nm), semi-minor axis (nm), orientation (deg N), and containment percent. This 
SPA represents a Bivariate Normal distribution for the target's location that has been 
obtained from an external source. [Ref. 3: pp. 12, 21] 

A. PATROLLING TARGET MODEL 

In the patrolling target model, the target is assumed to be randomly patrolling in the 
vicinity of the SPA. The target selects an initial course from a Uniform (0, 360) dis¬ 
tribution and an initial speed from a truncated triangular distribution. This triangular 
distribution takes as input the average target speed. The distribution ranges from half 
of the input speed to two times the input speed with the average speed being eight times 
as likely as cither of the endpoints. These distributions are shown in Figure 1 below. 



Figure 1. Patrolling Target Course and Speed Distributions 
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The target is assumed to continue on this course and speed for an amount of time 
determined from an exponential distribution with mean time between course changes 
being input by the user. At the completion of this leg in the target's track, an inde¬ 
pendent course, speed, and time on leg are selected from the same distributions and the 
process continues throughout the scenario. (Ref 3: pp. 10-11, 20| 

The values of the input parameters for a patrolling target were set as follows to ob¬ 
tain SALT output for testing: 

• Average target speed (knots) 4 

• Average time between course changes (hours) 3 

B. TRANSITING TARGET MODEL 

The transiting target motion model assumes that the target has a preferred course 
and speed of advance (SO.A) to travel. In this model, the user is required to input five 
parameters. The parameters are average target course (deg N). maximum variation in 
target course, average SOA (kts), maximum variation in SOA, and average time between 
course changes. For each leg. the target selects a course from a truncated triangular 
distribution ranging over the values average target course ± maximum variation in tar¬ 
get course with the average target course in the center being eight times as likely as ei¬ 
ther of the endpoints. The same type of distribution is used for selecting the target's 
speed of advance in the direction of the average target course on each leg. These dis¬ 
tributions are shown in Figure 2 below. 
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Figure 2. Transiting Target Course and Speed of Advance Distributions 


Each leg is assumed to be maintained for an exponential amount of time as dis¬ 
cussed in the Patrolling Target Model section above. (Ref 3: pp. 10-11, 20) 

I'he values of the input parameters for a transiting target weic set as follows to ob¬ 


tain SALT output for testing: 

• Average SOA (knots) 12 

• Maximum variation in SOA (knots) 4 

• Average target course (deg N) 045 

• Maximum variation in target course (deg) 30 

• Average time between course changes (hours) 6 

C. FLEEING DAI UM TARGET MODEL 


'I he model for fleeing datum target motion assumes that the taiget is clearing datum 
from an area where it suspects it has been detected. To this end, the target selects a 
course and speed from the specified distributions and maintains that course and speed 
throughout the encounter. The target couise is a draw from a Uniform (0, 360) distrib¬ 
ution. The target speed is selected from a truncated triangular disttibution with four 
input parameters. Ihe parameters to be input arc average speed, minimum speed, 
maximum speed, and average speed likelihood factor. The fust three parameters are 
self-explanatory, and the average speed likelihood factor determines how many times 
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more likelj the average speed is over the endpoint speeds. These distiibutions arc shown 
in Figure 3 on page 6. 



rigiire 3. Fleeing Datum Target Course and Speed Distributions 

Since the target maintains course and speed the entire time, there is no requirement for 
a mean time between course changes parameter. (Ref. 3: pp. 10-11, 20) 

T he values of the input parameters for a fleeing datum target wcic set as follows to 


obtain SALT output for testing: 

• .Average target speed (knots) 12 

• Maximum speed (knots) 15 

• Minimum speed (knots) 9 

• Average speed likelihood factor 8 


This completes the discussion of tlie three upcs of target motion models available 
for user selection in the SALT algotithm. T he next chapter will discuss the development 
and evaluation of sonobuoy patterns to be dcplovcd against each specific tv pc of target 


motion. 







III. SALT SONOBUOY PATTERNS 


The sonobuoy patterns recontniendcd by SALT for each type of target motion 
model are presented here in brief overview form. For each type of target motion, both 
the pattern development and evaluation of probability of detection of that pattern are 
discussed. The overall objective of the sonobuoy patterns developed by the initial search 
planning phase of SALT is to maximize the probability of detecting the target. 

A. PATROLLING TARGET PATTERNS 

A target which is patrolling in an area as discussed in the previous chapter will have 
an average velocity of zero and will tend to remain in the vicinity of the initial Search 
Probability Area (SPA). The optimal search density along each of the axes of the 
bivariate normal distribution would be parabolic. The optimization of search effort 
would require solving a very complicated nonlinear function and the density of search 
effort required would be very dilTicult to obtain with sonobuoys. In order to overcome 
the dilliculties associated with attempting to optimize the search effort, the SAL'f algo¬ 
rithm uses an approximation to the optimal search effort. The size of a rectangle which 
most closely appro.ximaics the bisariatc normal distribution of target location is first 
determined. The available search effort is then distributed uniformly over tliis rectangle. 
This process is described in more detail below. [Ref 2: p 36-37] 

As the search for the target progresses without any detections, the ' ‘liptical area of 
uncertaintv of the initial SPA continues to expand. To capture a static value for the si/c 
of the SP.A to use in detci mining and evaluating the sonobuoy dcplovn'icnt pattern, the 
initial SPA is first expanded to co: respond to the time which is half way througii the 
input search duration. This expansion gives larger values for a, and which repicscnt 
the average growth in uncertaintv in the targets position over the duration ol'tite search. 
(Ref. 2: p. 20-22] 

The bivariate normal distiihution dcteimincd by the SPA, expanded as described 
above, can bo expressed as the product of two independent univariate normal distrib¬ 
utions defined on the principal axes of the SP.A. Let a, and a. be llie standaui deviation 
of the target distiibution along the first and scco- 1 axes. The lengths of the sides of the 
rectangle p.irallel to tlicsc axes arc 2K<7i and 2 Krr 2 respectively for some constant K. 
These values determine a familv of rectangles with the parameter K detcimining tlie si/e 
of the rectangle as illustrated in 1 iguie 4 below. [Ref 2: pp. 37-38] 
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riguie 4. Family of Rectangle Plans 


To determine the probability of detection for a rectangular pattern of a given size, 
it is necessary to multiply the probability that the target is located in the lectangle by 
the conditional probability that the target is detected given that it is in the rectangle. 
1 his conditional probability is approximated by an exponential detection function which 
is based on the pessimistic assumptions of a random search (Ref 4: p. 28|. This idea 
of removing the conditioning is used by the initial search planning phase of SAL'f to 
deteimine that value of K for which the piobability of detection calculated by the fol¬ 
lowing c(|uation is maximized. 


iyK) = [dHK)--<I)(--K)]\l -e.xp(-f;-P^)J (3.1) 

4k (I|02 


where: 

n = the number of .sonobuoys 

W = (he sweep width of each sonobuoy in nautical miles 
V = average speed of target as input by the user 
r = search duration in hours 








Tilis maximum value is output by SALT as the probability of detection for the recom¬ 
mended sonobuoy pattern. [Ref 2: pp. 39-42] 

It should be noted here that the term nWVT in equation (3.1) can be thought of as 
the total area co\ercd by the pattern during time period T due to target motion through 
the sonobuoy field. We will later argue that a more accurate probability of detection 
can be calculated by accounting for the possibility of detecting the target at the begin¬ 
ning of the search. 

To develop this operationally feasible sonobuoy pattern, SALT uses the optimally 
sized rectangle calculated above and distributes the number of sonobuoys available over 
that rectangle to approximate a uniform distribution of search effort. To accomplish 
this distribution of effort, the sonolmoys are positioned in equally spaced rows witli the 
spacing between the buoys in a row being equal. If the number of sonobuoys is not di¬ 
visible bv the number of rows, remaining buoys arc allotted to rows beginning witli the 
center row and working outward. The sonobuoy pattern is output to the screen graph¬ 
ically along with the estimate of probability of detection. [Ref 2: pp. 43-50] 

An example of the graphical sonobuoy pattein output for a scenario involving 16 
sonobuoys deployed against a patrolling taigct is presented in Figure 5 on jiage 10. The 
inner ellipse represents the initial SPA and the outer ellipse represents the midscarch 
duration ellipse used to pcrfoim the above calculations. The optimally sized rectangle 
is also displayed in the pattern. 
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Figure 5. Recommended Pattern for 16 Sonobuoys vs. a, Patrolling Target 

B. TRANSITING TARGET PATT ERNS 

A target which is assumed to be transisting has a preferred direction of motion 
which can be exploited by the search planning algorithm by using barriers which "walk" 
against the target motion. This fact makes it possible for a line of sonobuoys to sweep 
out a large rectangular area. The size of the area swept depends on target speed in the 
direction of the barrier. The development and evaluation of sonobuoy patterns for de¬ 
ployment against a transiting target is similar to that discussed in the previous section 
for patrolling targets. The major dificrencc in the two developments is that the 
transiting target search plans are constructed in the relative motion space of the mean 
velocity of the target and then translated to geographic space. [Ref. 2: pp. 50-59] 
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An example of the graphical output of the recommended deployment pattern for 16 
sonobuoys against a transiting target with input parameters as discussed in Chapter II 
is shown in Figure 6 below. The ellipse displayed in the figure represents the initial 
SPA. The ellipse which is updated to the midsearch duration for use in performing cal¬ 
culations is not shown. 



Figuie 6. Recoiiiineiided Pattern for 16 Sonobuoys vs. a I'ransiting Target 

C. FLEEING DAFUM TARGET PATTERNS 

The development of a sonobuoy pattern to be dcplojcd against a ficcing datum 
target begins by updating the initial SPA to the mid-search time. This update is ac¬ 
complished using a discrete time and space Markov chain. As the search progresses, the 
distribution of target location spreads radially outward eventually becoming donut 
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shaped with the mean, or center, of the SPA becoming a low probability area. (Ref. 2: 

p. 62] 

In order to accommodate the distribution for a fleeing datum target, SALT uses a 
class of search plans referred to as "two rectangle plans". These plans are defined by an 
inner and an outer rectangle oriented concentrically such that the sides of the rectangles 
are parallel to the major and minor axes of the SPA. The total search eflbrt is spread 
uniformly over the annular region between these tw'o rectangles. The stated objective 
of SALT is to select the sizes of these rectangles so as to maximize the probability of 
detecting the target. (Ref. 2: p. 62] 

The probability of detection for these fleeing datum sonobuoy patterns is computed 
in a manner similar to the calculation discussed above for a patrolling target. The 
probability that the target is contained in the two rectangle annulus is multiplied by the 
conditional probability of detecting the target given that it is contained in this region. 
This conditional probability is determined by an exponential detection function as was 
done in the patrolling target case. [Ref. 2: p. 65] 

An example of this type of pattern for 16 sonobuoys as output by SALT is shown 
in Figure 7 on page 13. The sonobuoys in this pattern are located so far from the initial 
SPA due to the adjustment of an input parameter to the SALT scenario known as 
time-late which is the elapsed time between when the SPA was obtained and the time 
at which the search begins. In this scenario, time-late has been set to four hours. 

In performing the above calculations of probability of detection for the sonobuoy 
patterns to be depLyed against any of the three types of assumed target motion, SALT 
assumes that the sonobuoys are all dropped in place instantly at the beginning of the 
search. This assumption does not adversely aflect the results of either the patrolling 
target or ^he transiting target. In the case of the patrolling target, the target tends to 
remain in the search area for the duration of the search thus being exposed to the entire 
pattern as if it were laid instantaneously. The transiting target sonobuoy patterns are 
generally short, straight line patterns that do not take a long time to lay. This fact 
combined with the time required for the target to intercept the pattern tend to make this 
assumption plausible. 1 he sonobuoy patterns deployed against the fleeing datum target, 
however, are usually laid a large distance from the initial target SPA with large spacing 
between the sonobuoys. Due to the length of time required to lay these patterns and the 
higher speeds associated with fleeing datum targets, the assumption of having the 
sonobuoys in place initially is faulty. During the data collection runs of the fleeing da¬ 
tum target scenarios, the target was sometimes seen to be beyond the reach of the pat- 
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rigiire 7. Recommended Pattern for 16 Sonobuoys vs. a Fleeing Datum Target 


tern by tlie time the critical portion of the pattern was put in place. Since the goal of 
this thesis is only to test the validity of the probability of detection as output by SALT, 
this same assumption has been built into the simulations discussed in the following 
chapter. 
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IV. SIMULATION DOCUMENTATION 


Two examples of the FORTRAN simulation routines used to test the probability 
of detection of the sonobuoy patterns recommended by SALT as well as to determine 
if better patterns could be developed are presented in Appendix B. Seven FORTRAN 
programs are required to accomplish the desired testing. For each target motion model 
there are two separate programs. The first of these programs is written to be used wdth 
Definite Range Law sonobuoys, and the second program simulates sonobuoys with de¬ 
tection ranges that var\' according to a Lambda-Sigma Jump model. (SALT will handle 
either case.) The remaining program is a generic detection subroutine (DTCT) which is 
used by all of tlie simulations. Each of these programs is discussed below'. All of the 
simulation routines are developed using the SALT assumption that each sonobuoy is 
placed in its exact location at the beginning of the search. 

A scries of calibration runs has been executed for each of the target motion simu¬ 
lation programs discussed below to ensure proper operation. This testing was accom¬ 
plished by having the program write the endpoints of each track segment to an output 
file. This output was then graphed and compared to the respective target motion model. 
Each simulation show'ed target motion consistent with the assumptions for that type of 
target. 

A. PATROLLING TARGET SIMULATIO S 

The simulation routine for a patrolling targ with Definite Range Law sonobuoys 
is named "PATROL". The first step in the program is to read in the input parameters 
from a data file. The eight parameters needed by the simulation are as follows: Ox and 
Oy , the standard deviations of the normal distributions required to simulate the initial 
Search Probability Area (SPA) and its containment percent, search duration in hours, 
target's course change rate in inverse hours, detection range of the sonobuoys in nautical 
miles, an integer seed for use in the random number generator, average target speed in 
knots, and finally, the Cartesian coordinates of the location of each of the sonobuoys in 
the pattern. 

The target parameters arc simulated as follows. The initial target position is deter¬ 
mined by random draws from two independent normal distributions both with mean 
zero and standard deviations of Ox and Oy The target course, in radians, is selected from 
a Uniform (0, 2?:) distribution. The tiiangular distribution required for the target speed 
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is simulated by converting a Uniform (0, 1) random number into a speed using the in¬ 
verse CDF technique. The speeds generated by this technique are displayed as a 
histogram in Figure 8 on page 15. It can be seen from this figure that the speeds for 
patrolling targets in this simulation are in fact representative of a triangular distribution. 



'1 he time that the target will spend on the chosen course at tlie chosen speed is deter¬ 
mined by a random draw from an e.xponcntial distribution with mean 
(course change rate) ’ . These three parameters, course, speed, and time, determine a 
line segment uhich repiesents a leg of the target track. The final target position for this 
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leg of the track is calculated using the initial position, course, speed, and time on leg 
determined as discussed above. 

The coordinates of a sonobuoy in the pattern, the detection range of that sonobuoy, 
and the endpoints of the current leg of the target track are used in a call to the DTCT 
subroutine to determine whether or not a detection has occurred. This process is re¬ 
peated until all sonobuoys in the pattern under investigation have been checked for 
possible detections. If a detection has occurred, the time of the detection is calculated, 
and the number of detections is incremented by one. The detection time of that 
sonobuoy which detects the target first in the time sequence is written to an output file 
for use in further analysis, and this run of the simulation is terminated. If no detections 
take place on this leg, the final position of this leg is reset to be the initial position on 
the next leg and the time that the search has been underway is incremented. If this time 
is less than the search duration input by the user, a new leg and final position are se¬ 
lected as above and the process is repeated. 

The program which simulates a patrolling target in a sonobuoy field in which de¬ 
tection range is fluctuating according to a Lambda-Sigma Jump model is entitled 
"LSPAT". The input parameters for this program are the same as those discussed above 
with one exception. Instead of a constant value of sonobuoy detection range being in¬ 
put, a Figure of Merit (FO.VI) is input by the user in dB. 

The target parameters, course, speed, and time on leg, are calculated as described 
above. However, a new environmental parameter is introduced into this model. This 
parameter is the acoustic fluctuation in dB which is represented by a random draw from 
a Normal distribution with mean zero and standard deviation, sigma = 8 dB. This 
acoustic fluctuation remains constant for a period of time determined from an expo¬ 
nential distribution with mean, lambda = 1 hour. 

1 he acoustic fluctuations simulated as discussed in the previous paragraph are con¬ 
sidered to be ambient acoustic fluctuations which apply uniformly to all sonobuoys in 
the pattern. The latest acoustic fluctuation is added to the input FOM to determine the 
current FOM. The current FO.M is used to determine the detection range of the 
sonobuoys by linear interpolation between the two nearest points on the transmission 
loss curve which is input as a separate data file. The detection range calculated in this 
way applies to all of the sonobuoys in the pattern. 

Calls to the DTCT subroutine for each sonobuoy to dctc.''mine detection,non- 
detection are treated as discussed abose. I he only diflerence is that the current segment 
of the target track has its final position either at the end of the acoustic fluctuation time 
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or at the end of the time on the current leg whichever occurs next in the time sequence. 
If the completion of acoustic fluctuation time occurs first, then time on leg is decre¬ 
mented by the acoustic fluctuation time and total elapsed search time is incremented by 
the same amount of time. A new acoustic fluctuation and time for the fluctuation are 
then drawn . If the completion of the time on leg occurs first, then the acoustic fluctu¬ 
ation time is decremented by leg time and total time is incremented by it. New values 
of course, speed, and time on leg are then computed. 

If detection occurs, the number of detections is incremented by one, the time of the 
detection is written to an output file, and this run of the simulation is finished. If no 
detections have occurred and the total elapsed search time is less than the search dura¬ 
tion input by the user, the process continues. 

B. TR.4NSmNG TARGET SIMULATION 

The program which simulates a transiting target versus a field of Definite Range 
Law sonobuoys functions in a manner similar to that described in the above section on 
Patrolling Target .Simulations. The input parameters for this program are the same as 
those discussed above except that the minimum and maximum course are input in 
.adians and the minimum and maximum speed arc input in knots. 

The transiting target parameters are simulated as follows. The initial position is 
determined the same way as described under the Patrolling Target Simulations section 
above. Both target course and speed of advance (SO.A) are simulated from triangular 
distributions using the inverse CDF technique. These values of course and SOA are then 
converted to target speed by dividing the SOA by the cosine of the angle between aver¬ 
age target course and the target course determined for this track segment. These trian¬ 
gular distributions have been verified by histogram as in the patrolling target case, but 
the figure is not reproduced here. The remainder of the simulation proceeds exactly as 
described above for patrolling targets. 

The program which models the Lambda-Sigma Jump sonobuoys generalizes from 
the transiting target Definite Range Law sonobuoy program in the same manner as for 
patrolling targets. 

C. FLEEING DATUM TARGET SIMULATION 

The next program to be discussed is the one used to simulate fleeing datum targets 
and Definite Range Law sonobuoys. The input parameters are the same as those dis¬ 
cussed in the previous sections. The target course is drawn from a Uniform (0, In) dis¬ 
tribution and target speed is drawn from a triangular distribution determined by the 
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input parameters of minimum and maximum target speed. The course and speed need 
only be drawn one time for each simulation run since the target is assumed to maintain 
that course and speed throughout the engagement. The initial position is aged by the 
time-late value of four hours using the selected course and speed. The deteccion/non- 
detection event is then determined as described in the Patrolling Target Simulations 
section. This simulation routine generalizes to the Lambda-Sigma Jump model for 
sonobuoy detection range in the same manner as previously discussed. 

D. DETECTION SUBROUTINE 

The detection subroutine entitled "DTCT" is used by all six of the main programs 
to determine whether or not a target detection occurs. If a detection does take place, 
the subroutine returns the distance from the beginning of the current segment of the 
target track to the point where the earliest detection occurs. 

The inputs to the DTCT subroutine aie the detection range of the sonobuoys. the 
cartesian coordinates of the endpoints of the current track segment, and the location of 
the current sonobuoy being tested for detection criterion. These seven parameters are 
passed into the subroutine at each call to it. 

The first step of this subroutine is to calculate the perpendicular distance between 
the line determined by the target's track segment and the location of the sonobuoy. If 
this value is greater than the detection range of the sonobuoy, no detection takes place. 
However, if this perpendicular distance is less than the detection range, there arc three 
cases to consider to determine wnethcr -vr not detection occurs. If the track segment 
contains the point on the line that is closest to the buoy center, then a detection occurs. 
If the track segment docs not contain that point, but either endpoint of the track seg¬ 
ment falls within the range of detection of the buoy a detection will also take place. If 
neither of these conditions holds, then there is no detection for this track segment- 
sonobuoy combination. Each of the.se four cases is illustrated in Figure 9 on page 19. 

In order to ensure the proper operation of this subroutine, examples of the four 
cases described above were developed using simple graphical techniques. These exam¬ 
ples were then input into the subroutine for testing. The outputs from the subroutine 
confirm that it performs its function properly as described here. 
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Tigure 9. 





Soiiobuoy Detection/Noii-detectlon Possibilities 
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V. SIMULATION RESULTS OF SALT SONOBUOY PATTERNS 

The sonobuoy patterns generated by the SALT algorithm for each specific type of 
target motion were tested using the respective simulation routine discussed in the previ¬ 
ous chapter to determine if the probability of detection for that pattern as output by 
SALT is accurate. The simulation results are output as a number of detections which 
can be easily converted to a probability of detection if the number of detections is di¬ 
vided by the number of runs of the simulation which is 10,000 runs in each case. Each 
proportion, p, which represents the probability of detection can be expanded to a two- 
sided (I — a)% Confidence Interval by use of the following formula: 


P 


+ 7,._a_ 
2 


Pji-P) 

n 


(5.1) 


where Z,_^ represents the proper percentage point of a Standard Normal cumulative 
distribution. 

The SALT algorithm was run three times for each type of target motion using 
sonobuoy s that have a Definite Range Law detection function. The parameters within 
the SALT program which determine the acoustic fluctuations were set to their minimum 
values to accomplish this (see Appendix A). For each run of a specific target motion 
type, all inputs remain the same with the exception of the number of sonobuoys. The 
number of sonobuoys used for each run are 8, 16, and 24. These numbers of sonobuoys 
allowed for testing over a range of values of probability of detection. The same nine 
runs of the S.4LT algorithm were repeated to obtain data for sonobuoys which behave 
according to a Lambda-Sigma Jump model. The parameters for this phase of testing 
are an ambient acoustic fluctuation rate of one per hour and standard deviation of am¬ 
bient acoustic fluctuations of eight dB. The probability of detection as output by SALT 
will now be compared to the simulation results for each type of target motion. 

A. PATROLLING TARGET 

The probability of detection for the Definite Range Law sonobuoy patterns recom¬ 
mended by S.ALT for a patrolling target are displayed in Table 1 on page 21. Both the 
S;\LT calculated probability of detection and the simulation routine output probability 
of detection are listed for the number of sonobuoys used in each pattern. A 95Vo Con- 
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fidence Interval, calculated as described above, is also displayed for each simulation 
probability of detection. 


Table 1. PATROLLING TARGET COMPARISONS (DEFINITE RANGE) 


NUMBER 

SALT 

SIMULATION 

95% CONFIDENCE 

OF BUOYS 

Pd 

P D 

INTERVAL 

8 

0.49 

0.8848 

(0.8785. 0.8911) 

16 

0.67 

0.9402 

(0.9378, 0.9426) 

24 

0.76 

0.9773 

(0.9758, 0.9788) 


It can be seen from this table that the SALT algorithm tends to understate the 
probability of detection for a target in this case. The low probabilities of detection 
output by S.ALT have the potential to be very inisleading. For example, the SALT re¬ 
commended pattern for eight sonobuoys lists as its probability of detection a value of 
0.49. This low probability of detection may encourage a P-3C crew to opt instead for a 
sonobuov pattern with more buo\s which would be a waste of resources considering the 
actual probability of detection for the target using an eight buoy pattern is 0.88. 

The Technical Documentation of SALT claims that the exponential detection func¬ 
tion used in equation (3.1) on page 9 gives a lower bound on the probability of detection 
that can be achieved by a systematic uniform search of a region (Ref 2; p. 26]. Due to 
this fact, one would expect to get resulting probabilities of detection that arc a little 
higher than the SALT output values. However, the disparity in the values given here is 
too high to be so easily e.xplained away. In order to gain insight into why the simulation 
probability of detection is so different from the SALT probability of detection, the de¬ 
tection times as output by the simulation program are graphed into a Cumulative Dis¬ 
tribution Function in Figure 10 on page 22. 
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DEJECTION TIME (HOURS) 


Figure 10. Cumulative Distribution of Detection limes 

It can be seen from this CDF that a large percentage of the targets whicli get de¬ 
tected arc initial detections. This is not too surprising considering the tight packing of 
tlic sonobuoys into the containment ellipse. The sonobuoy detection circles are 
shown plotted over this 86®o containment ellipse in Figure 11 on page 23. The ellipse 
sliown in this figure has a semi-major axis of 60 nautical miles and a semi-minor axis of 
30 nautical miles. Ihe detection range circles are 18.2 nautical miles in radius and are 
placed in the locations specified by SALT. 


22 




Figure 11. Detection Ranges and 86% Containment Ellipse for an 8 Sonobuoy 
Pattern 

To account for this initial probability of detection in the patrolling target sonobuoy 
patterns, equation (3.1) on page 9 can be modified slightly to read as follows; 

Pd(K) = [<b(K) - (I)( -K)]'[l - exp( - )] (5.2) 

4K a,02 

where: 

R = the detection range of the sonobuoys in nautical miles 

and all other variables arc defined as they arc in equation (3.1). 

fhe above modification to the calculation of the probability of detection takes into 
account the detection area of each sonobuoy at the initiation of search, 'fhe value of 
the product nWVf is clearly time dependent, 'fhis means that in equation (3.1) there 
is zero probability of an initial detection. The detection times shown in Figure 10 on 
page 22 clearly show the need to account for this initial piobability of detection. 







Using equation (5.2) to calculate the probability of detection for the three Definite 
Range Law sonobuoy patterns specified in Table 1 on page 21 results in probabilities 
of detection of 0.74, 0.89, 0.93 for 8, 16, and 24 sonobuoy patterns respectively. While 
these values do not exactly reflect the results of the simulation runs, they do provide the 
decision maker with a much more reahstic assessment of the probability of detection for 
the recommended pattern. 

The Lambda-Sigma Jump model comparisons for the similar patterns are given in 
Table 2. The results for the 8 and 16 sonobuoy patterns show the same phenomenon 
discussed above for Definite Range Law sonobuoys. In these cases, equation (5.2) again 
provides a more realistic value of probability of detection. 


Table 2. PATROLLING TARGET COMPARISONS (LAMBDA-SIGMA) 


NUMBER 

SALT 

SIMULATION 

95% CONFIDENCE 

OF BUOYS 

Pd 

Pd 

INTERVAL 

8 

0..59 

0.8664 

(0.8630, 0.8698) 

16 

0.76 

0.9028 

(0.8994. 0.9062) 

24 1 

0.84 

0.8120 

(0.8081.0.8159) 


The SALT output for 24 sonobuoys in this scenario was a pattern which was ro- 
ticeably different from the five other sonobuoy patterns recommended for deployment 
against a patrolling target. The S.ALT generated pattern is a BRUSHTAC pattern where 
all of the buoys are placed in a narrowly spaced line along the major axis of the initial 
SPA. In a telephone conversation with Lawrence D. S^one, the president of Metron, 
Inc., it has been determined that an error in the SALT code resulted in an optimal target 
depth being selected which would give the sonobuoys a convergence zone detection. 
This selection was made by the algorithm in spite of the fact that the probability of the 
target being at that depth was input as zero. Due to this unfortunate fact, the results 
displayed in Table 2 for a 24 sonobuoy pattern have no meaning. 

B. TRANSITING TARGET 

Ihe comparison of the SALT output probability of detection for transiting target 
sonobuoy patterns to the probability of detection for the same patterns as output by the 
simulation routine is presented in Table 3 on page 25. 


24 

























Table 3. TRANSITING TARGET COMPARISONS (DEFINITE RANGE) 


NUMBER 

OF BUOYS 

SALT 

Pd 

SIMULATION 

P D 

95% CONFIDENCE 
INTERVAL 

8 

0.72 

0.8857 

(0.8825, 0.8889) 

16 

0.81 

0.9085 

(0.9056, 0.9114) 

24 

0.82 

0.9156 

(0.9128,0.9184) 


The probabilities of detection displayed above do not exhibit the same magnitude 
of difference that those of the patrolling target exhibit in the previous section as seen in 
Table 1 on page 21. However, there is still a slight difference in these probabilities of 
detection that requires consideration. The CDF plot for transiting target detection times 
reveals that only about 35% of the detections obtained are due to initial detection. The 
reason for this smaller number of initial detections is that the transiting target search 
plans are constructed in the relative motion space of the mean velocity of the target as 
discussed in the Transiting Target Patterns section of Chapter III. The length of the side 
of the rectangular area swept out which is parallel to the mean velocity vector is deter¬ 
mined by the following equation; [Ref 2: p. 52) 

Length = VT 2R (5.3) 


where: 

V = average SOA of the target in knots 

T = search duration in hours 

R = the detection range of the sonobuoys in nautical miles 

When this relative motion space is translated into geographic space, only the 2R 
portion of the length as determined by equation (5.3) has the possibility for initial de¬ 
tections. In order to account for these initial detections in the calculation of probability 
of detection, this fraction of the total length should be used. 

Due to the same depth selection error discussed at the end of the previous section, 
the SALT output for transiting targets using the Lambda-Sigma Jump model is only 
valid for the eight sonobuoy case. The SALT output probability of detection for this 
scenario is O.Sl while the simulation output for the identical pattern is 0.7996. The 95” o 



















Confidence Interval about this value is (0.7956, 0.8036). In this more realistic case, the 
initial probability of detection does not seem to be an important factor. 

C. FLEEING DATUM TARGET 

The results of the simulation analysis of Definite Range Law sonobuoy patterns 
deployed against a fleeing datum target are shown in Table 4. 


Table 4. FLEEING DATUM TARGET COMPARISONS (DEFINITE RANGE) 


NUMBER 

OF BUOYS 

SALT 

Pd 

SIMULATION 

Pd 

95% CONFIDENCE 
INTERVAL 

8 

0.47 

0.4829 

(0.4779, 0.4879) 

16 

0.69 

0.9579 

(0.9559, 0.9599) 

24 

0.81 

0.9985 

(0.9981,0.9989) 


The probability of detection as output by both SALT and the simulation of this 
scenario appear to be in close agreement for the eight sonobuoy case while the other two 
patterns show a highly significant difference. This difference cannot be explained by 
initial detections as is done in the cases of patrolling and transiting targets above. Less 
than 10°/o of the detections were achieved initially when these scenarios were run using 
the simtilation routines. An alternate reason for this phenomenon can be found by de¬ 
termining the fraction of the perimeter of these patterns that is co\ered by sonobuoy 
detection radii. In the eight sonobuoy pattern, 48“/o of this perimeter is covered by such 
detection ranges. Similar calculations for the 16 sonobuoy pattern reveal that 94" o of 
the perimeter is so covered, and the percentage of coverage for the 24 sonobuoy case is 
100"'o. This fact accurately reflects that when the sonobuoys are assumed to be placed 
in the water immediately upon commencement of search, it is highly unlikely that a 
target will avoid having to pass through the annular region where the sonobuoys are 
placed. 

As discussed in the Fleeing Datum Target Patterns section of Chapter III, the 
probability of detection in these scenarios is calculated within the SALT algorithm by 
multiphing the probability of detecting the target given that it is contained in the two 
rectangle annular region by the probability that the target is contained in this region. 
By performing this simplified calculation, SAL f fails to take into account the fact that 
a target which is inside of the inner rectangle at the commencement of search must pass 
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through this annular region to evade detection. The SALT algorithm must be modified 
in some way to account for this possibility in order to more realistically evaluate the 
probability of detection for fleeing datum targets. 

The results for the same target scenarios versus Lambda-Sigma Jump model 
sonobuoys are presented in Table 5. These results show a disparity similar to that dis¬ 
cussed above for Definite Range Law sonobuoys. Thus the corrective action specified 
above is valid for both types of sonobuoy detection models. 


Table 5. FLEEP'G DATUM TARGET COMPARISONS (LAMBDA-SIGMA) 


NUMBER 

OF BUOYS 

SALT 

SIMULATION 

Pd 

95% CONFIDENCE 
INTERVAL 

8 

0.59 

0.8623 

(0.8589. 0.8657) 

16 

0.80 

0.9138 

(0.9110, 0.9166) 

24 

0.88 

0.9413 

(0.9389, 0.9437) 


It can be seen from these analyses that in all cases of target motion models the 
SALT code requires modification to be able to more accurately reflect the probability 
of detection so that sound decisions can be made by the user of this tactical decision aid. 
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VI. ALTERNATIVE SONOBUOY PATTERN DEVELOPMENT 


The second phase of this research effort involves attempting to determine whether 
or not the sonobuoy patterns output by the SALT algorithm can be improved. The 
measure of effectiveness of a sonobuoy pattern will be the probability of detection as 
output by the same simulation routines used to verify the probability of detection. The 
improvements considered and the results are presented for each target motion model in 
the following sections. 

A. PATROLLING TARGET MODEL 

As was shown previously in Figure 11 on page 23, the pattern for eight r.onobuoys 
deployed against a patrolling target consists of excessive overlap of the buoy detection 
regions. The sonobuoy spacing in this pattern is 22.3 nautical miles between the two 
rows and 21.2 nautical miles between the buoys within a row. The probability of de¬ 
tection of this pattern is 0.88 as determined from the simulation program. 

The SALT assumption for deploying a sonobuoy pattern against a target of this 
type is that the search effort should be distributed uniformly over the region in which the 
target is most likely to be located. Under this assumption it can plainly be seen that the 
overlap of the buoy detection ranges is excessive. By increasing the spacing between the 
buoys in this pattern, some gains in probability of detection can be realized. 

The best spacing to use between the sonobuoys in the pattern with eight Definite 
Range Law sonobuoys has been determined by trail and error. After many simulation 
runs on patterns with various spacing, it has been determined that the best effort can 
be achieved by making the spacing between rows 31.2 nautical miles, and setting the 
within row spacing to 36.0 nautical miles. The probability of detection for this pattern 
was output by the simulation as 0.96. Further spreading of the sonobuoys results in too 
much effort being allocated to outlying areas. The probability that the target is located 
in these areas is very low. 

A graphical representation of the detection ranges of the sonobuoys superimposed 
on an 86” o probability of containment ellipse is shown in Figure 12 on page 29. The 
ellipse shown in the figure has a semi-major axis of 60 nautical miles and a semi-minor 
axis of 30 nautical miles. The detection range circles arc 1S.2 nautical miles in radius, 
and overlap only a slight amount. 
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Figure 12. Detection Ranges and 86% Containment Ellipse for an 8 Sonobuoy 
Pattern nith 36 Nautical Mile Spacing 

Similar tests were run using the Lambda-Sigma Jump model sonobuoys, but altering 
the spacing does not alFect the results as much since there are changes in sonobuoy de¬ 
tection range caused by acoustic fluctuations. The SALT output spacing for the eight 
sonobuoy Lambda-Sigma Jump pattern is 24.7 nautical mile.s between rows and 23.5 
nautical miles between buoys within each row. Altering the spacing between the 
sonobuoys in this scenario provides no significant increases in the probability of de¬ 
tection. 

Another type of pattern that was tested using Lambda-Sigma Jump model 
sonobuoys against a patrolling target is a 2 x 4 x 2 pattern. This pattern consists of 4 
sonobuoys placed along the major axis of the 86" o containment ellipse with 2 sonobuoys 
on each side of this line located near the minor axis. No notable increase in probability 
of detection is realized using this pattern. 




Based on the results of the testing described above, it seems that the sonobuoy pat¬ 
terns recommended by SALT for use in the more realistic environmental setting modeled 
by the Lambda-Sigma Jump process are consistent with the stated objectives of the 
SALT algorithm. 

B. TRANSITING TARGET MODEL 

In attempting to improve the sonobuoy patterns recommended by SALT for use 
against a transiting target, only the more realistic case of the Lambda-Sigma Jump 
model sonobuoys is considered. The SALT pattern for this type of target motion model 
is a straight line of buoys placed perpendicular to the mean target .:ourse. In looking 
at this pattern of eight sonobuoys graphed along with the 86% containment ellipse, it 
appears that the line of sonobuoys is located too close to the center of the SPA. Moving 
the line of sonobuo> s away from the center in the direction of the mean target course 
results in the probability of detection going down. This results from the fact that al¬ 
though fewer targets will be located on the wrong side of this line, more targets will be 
located further from the pattern and will not reach it prior to termination of the search. 
The next aspect of the pattern to be considered is the spacing beween the sonobuoys 
in the line. The spacing between the sonobuoys in the SALT generated pattern is only 
12.8 nautical miles, even though the average detection range is 26.3 nautical miles. The 
sonobuoy spacing was increased to 20 nautical miles. The simulation run for this new 
pattern showed no notable increase in probability of detection. This result is not diffi¬ 
cult to understand in light of the fact that as the spacing between the sonobuovs is in¬ 
creased, the buo\s near the ends of the line are being moved to areas with extremely low 
probability of achieving a detection. 

The next logical step in the development of alternative patterns is to depart from the 
requirement to lay the sonobuoys in a straight line. By increasing the spacing between 
the sonobuojs to overcome the excessive overlap and then bringing the endpoints of the 
line closer to the SPA, it seems that the probability of detection should increase. Two 
patterns were tested using this theory. 

The first pattern tested in which the sonobuoys are not placed in a straight line is a 
semi-circular pattern with spacing between the sonobuoys of 20 nautical miles. This 
pattern covers 90 degrees of arc on each side of the mean target course line. The pattern 
is positioned so that all sonobuoys are located at a distance midway between the lengths 
of the semi-major and semi-minor axes from the center of the SPA. 
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The second pattern tested is a chevron pattern. The sonobuoys in this pattern are 
divided equally between the two sides of the chevron, and the wedge angle of the chevron 
is set at 90 degrees. The point of the chevron is located along the mean target course 
line at a distance half way between the lengths of the semi-major and semi-minor axes. 
The spacing between the sonobuoys w'as again set to 20 nautical miles. 

Neither of these patterns results in a significant increase in the probability of de¬ 
tection for the transiting target modeled here. Based on the test results described here, 
the sonobuoy patterns generated by SALT for deployment against transiting targets 
appear to adequately meet the objectives of SALT. 

C. FLEEING DATUM TARGET MODEL 

The first sonobuoy pattern output by SALT that is considered for improvement is 
the pattern with eight Definite Range Law' sonobuoys. The analysis of the detection 
times output by the simulation routine for this pattern show's only 42 initial detections 
out of the 4829 detections obtained. This result leads to the suspicion that this pattern 
is located too far from the initial SPA, Each side of the square pattern generated by 
SALT is located 88 nautical miles from the center of the SPA. Several runs of the sim¬ 
ulation have been conducted moving the sonobuoy pattern nearer the SPA each time 
and scaling it accordingly. These runs show that a pattern in which the sides are located 
50 nautical miles from the center increases the probability of detection to 0.78. 

The square shape of the sonobuoy patterns output by SALT for attempting to de¬ 
tect a fleeing datum target naturally suggests that a better use of available search effort 
would be to arrange the sonobuoys in a circle so that the buoys would be equidistant 
from the SPA center. Such a pattern was tested using the eight sonobuoy Lambda- 
Sigma Jump model scenario. The sonobuoys in this pattern were placed 82.5 nautical 
miles from the center of the SPA and arranged so as to uniformly cover the perimeter 
of the circle. This distance from the center was arrived at by averaging the semi-major 
and semi-minor axes of the initial SPA and then adding a distance to account for the 
time-late of the beginning of the search which is four hours. In this way, all sonobuoys 
in the pattern are located closer to the center than in the SALT pattern. Each edge of 
the square SALT pattern is located at 84.6 nautical miles from the center. 

A similar circular pattern has also been tested in which the sonobuoys are placed 
60 nautical miles from the center. This distance is based on the improved pattern for 
Definite Range Law' sonobuoys discus.sed above. This pattern only increases the prob¬ 
ability of detection by 0.03 which is not significant. 
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The above results of attempting to improve the sonobuoy patterns recommended 
by SALT again show that the Definite Range Law scenario patterns can be improved 
greatly while the Lambda-Sigma Jump scenario patterns do appear to be consistent with 
the SALT objectives. These results should not be weighted too heavily as they are ob¬ 
tained by remaining within the constraints of the SALT assumptions, most notably that 
all of the sonobuoys in a pattern deployed against a fleeing datum target are put in place 
at the commencement of the search. This assumption has been shown to be faulty at 
the end of Chapter 111. 
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VII. CONCLUSIONS 


It has been established at the beginning of this thesis that the goal of the research 
effort is to investigate the initial search planning phase of the Search and Localization 
Tactical Decision Aid. This investigation consists of first attempting to ensure that the 
probability of detection of various sonobuoy patterns recommended by SALT is correct 
under the stated assumptions and subsequently of testing alternative sonobuoy patterns 
which could lead to a better patterns than those recommended by SALT. 

The probability of detection output by SALT for its own patterns is significantly 
understated for both the patrolling target model and the fleeing datum target model. 
Possible corrections to these errors have been presented here. The probability of de¬ 
tection for the transiting target model appears to be correct especially in the more real¬ 
istic case of the Lambda-Sigma Jump model for sonobuoy detection ranges. 

The sonobuoy patterns generated by SALT for Definite Range Law sonobuoys de¬ 
ployed against patrolling and fleeing datum targets have been shown to be inadequate 
for accomplishing the goal of maximizing probability of detection. The recommended 
patterns for all three motion models when the acoustic fluctuations which influence 
sonobuoy detection range are modeled by a Lambda-Sigma Jump process arc consistent 
with the stated objectives of S.ALT. This fact is especially encouraging since the patterns 
output by the SALT algorithm are operationally easy to deploy from a P-3C aircraft. 
One must be reminded here that the testing performed on the sonobuoy patterns output 
by SALT has been conducted within the constraints of the assumptions made in devel¬ 
oping the S.ALT algorithm. 

One aspect of the fleeing datum target scenario which warrants further investigation 
is the assumption that all sonobuoys are put in place simultaneously at the commence¬ 
ment of the search. Some method of calculating the probability of detection for this type 
of target needs to be developed which allows for the fact that a target may be outside 
of the sonobuoy pattern by the time that the critical portion of these large patterns is 
put in place. A more useful alternative to altering the calculation of the probability of 
detection would be to develop a means of accounting for the required flight time in re¬ 
commending the sonobuoy pattern. This change in the algorithm would open up the 
possibility of deploying patterns which spiral out from the center as time elapses and the 
target continues its assumed motion away from the SPA center. |Rof 4; pp. 1 l.^'-l 18j 
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The SALT algorithm is rated as effective in its initial search planning phase for pa¬ 
trol'ing targets with only minor modifications necessary to improve the calculations of 
probability of detection. The SALT recommendations from the initial search planning 
phase for transiting targets also appear to be realistic. Further development is required 
in the case of fleeing datum targets as discussed above in order to make SALT ready for 
operational testing. 
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APPENDIX A. SALT INPUTS FOR TEST CASES 


The purpose of this appendix is to provide a list of the inputs to the SALT algorithm 
which were used to generate the sonobuoy patterns required for the testing performed 
here. Only those parameters which are relavent to the scenarios are presented. The re¬ 
maining inputs are left at their default values. The target motion model parameters are 
not presented here since they are given in Chapter II. 

The following list shows input values used in the Definite Range Law sonobuoy 
cases for patrolling and transiting targets. 


Setup 


Buoy acoustic fluctuation rate (per hr) 

0.01 

Standard deviation of buoy acoustic fluctuations (dB) 

0.01 

Ambient acoustic fluctuation rate (per hour) 

0.01 

Standard deviation of ambient acoustic fluctuations (dB) 

0.01 

Search duration (hrs) 

4.0 

Environmental Data 


Ocean 

Pacific 

Ambient noise level (dB) 

75 

Target Parameters 


Depth 1 (ft) 

90 

Probability of depth 1 

1.0 

Depth 2 (ft) 

300 

Probability of depth 2 

0.0 

Frequency 1 (hz) 

50 

Source level 1 (dB) 

160 

Frequency 2(hz) 

300 

Source level 2 (dB) 

75 
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SPA Data 

Latitude 50 00 N 

Longitude 020 00 W 

Senii-major axis (nm) 60.0 

Semi-minor axis (nm) 30.0 

Orientation (deg N) 90.0 

Containment percent 86.0 

The input parameters which are changed to generate output patterns for fleeing da¬ 
tum targets is presented in the next list. 

Setup 

Search duration (hrs) 8.0 

SPA Data 

Semi-major axis (nm) 30.0 

Semi-minor axis (nm) 15.0 

The two input parameters that are changed to simulate sonobuoys with detection 
ranges which vaiy according to the Lambda-Sigma Jump model are as follows: 

Setup 

Ambient acoustic fluctuation rate (per hour) 1.0 

Standard deviation of ambient acoustic fluctuations (dB) 8.0 
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APPENDIX B. SIMULATION PROGRAMS 


* This appendix presents listings of examples of the FORTRAN computer code used 

in the simulation of the sonobuoy patterns deployed against the three target motion 
» models of SALT. 

The two programs used to simulate patrolling target motion models are shown here. 
The first program shows the FORTRAN code w'hich tests the SALT output patterns 
generated using Definite Range Law sonobuoys. The second program provides a model 
for sonobuoys which behave according to a Lambda-Sigma Jump process. The pro¬ 
grams used to simulate the transiting target and the fleeing datum target are similar ex¬ 
cept in the required target input parameters. 

Following the patrolling target simulation listings is the listing of the detection 
subroutine. This subroutine is used by all six of the simulation programs to determine 
whether or not the sonobuoy detection criterion is met. 
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1 PROGRAM PATROL 

2 REAL SIGMAX, SIGMAY, SCHTM, LMBDA, DETRNG 

3 INTEGER SEED, DET, NUMDET, I 

4 REAL AVGSP, SPEED, U, COURSE, XI, Yl, X2, Y2, XC, YG 

5 REAL T, LEGTM, TOTTIM, DSTDET, DT, DETTIM 

6 CALL EXCMSC'FILEDEF 01 DISK PATPRM DATA A') 

7 CALL EXCMSC'FILEDEF 02 DISK PATBYS DATA A') 




8 


READCOl,*) SIGMAX, SIGMAY, SCHTM, LMBDA, DETRNG, SEED, AVGSP 



9 


NUMDET = 0 



10 


DO 40 I = 1, 10000 , 


1 

11 


CALL IjqORMCSEED, XI, 1, 1, 0) 


1 

12 


XI = XI * SIGMAX 


1 

13 


CALL LNORMCSEED, Yl, 1, 1, 0) 


1 

14 


Yl = Yl * SIGMAY 


1 

15 


TOTTIM = 0. 0 


1 

16 


DETTIM = 10000.0 


1 

17 

10 

CONTINUE 


1 

18 


CALL LRNDCSEED, COURSE, 1, 1, 0) 


1 

19 


COURSE = COURSE * 6.283185 


1 

20 


CALL LRNDCSEED, U, 1, 1, 0) 


1 

21 


IF CU .LE. 1. 0/3.0) THEN 

1 

1 

22 


SPEED = AVGSP*C24 + SQRTC3024*U + 16))/56 

1 

1 

23 


ELSE 

1 

1 

24 


SPEED = AVGSP*C60 - SQRTC1528 - 1512*U))/28 

1 

1 

25 r 


END IF 


1 

26 


CALL LEXPNCSEED, T, 1, 1, 0) 


1 

27 


LEGTM = T/LMBDA 


1 

28 


X2 = XI + SPEED * LEGTM * SINCCOURSE) 


1 

29 


Y2 = Yl + SPEED * LEGTM * COSCCOURSE) 


1 

30 

20 

CONTINUE 


1 

31 


READC02,*) XC, YC 


1 

32 


IF CXC .GT. 999.0) GO TO 30 


1 

33 


CALL DTCTCDETRNG, XI, Yl, X2, Y2, XC, YC, DET, DSTDET) 


1 

34 


IF CDET .EQ. 1) THEN t 

1 

1 

35 


DT = TOTTIM + DSTDET/SPEED 

1 

1 

36 


IF CDT .LT. DETTIM) THEN 

2 

1 

37 


DETTIM = DT 

2 

1 

38 


END IF 

1 

1 

39 


END IF 


1 

40 


GO TO 20 


1 

41 

30 

CONTINUE 


1 

42 


REWIND 02 


1 

43 


XI = X2 


1 

44 


Yl = Y2 


1 

45 


TOTTIM = TOTTIM + LEGTM 


1 

46 


IF CTOTTIM .LT. SCHTM .AND. DETTIM . GT. 9999.0) GOTO 10 


1 

47 


IF CDETTIM .LT. SCHTM) THEN 

1 

1 

48 


NUMDET = NUMDET + 1 

1 

1 

49 


WRITECll,*) DETTIM 

1 

1 

50 


END IF 


1 

51 

40 

CONTINUE 



52 


PRINT*, 'NUMBER OF DETECTIONS =', NUMDET 



53 


STOP 



54 


END 


3R 









1 


PROGRAM LSPAT 




2 


REAL SIGMAX, SIGMAY, SCHTM, LMBDA, DETRNG, FOM 



3 


REAL ACFL, AFTIM, R, Rl, TL, TLl 




4 


INTEGER SEED, DET, NUMDET, I 




5 


REAL AVGSP, SPEED, U, COURSE, XI, Yl, X2, Y2, 

XC, YC 



6 


REAL T, LEGTM, TOTTIM, DSTDET, DT, DETTIM 




7 


CALL EXCMSCFILEDEF 01 DISK PATPRM DATA A') 
CALL EXCMSCFILEDEF 02 DISK PATBYS DATA A’) 




8 





9 


CALL EXCMSCFILEDEF 03 DISK TLCRV DATA A') 




10 


READCOl,*) SIGMAX, SIGMAY, SCHTM, LMBDA, FOM, 

SEED, AVGSP 



11 


NUMDET = 0 




12 


DO 110 1=1, 10000 



1 

13 


CALL LNORM(SEED, XI, 1, 1, 0) 



1 

14 


XI = XI * SIGMAX 



1 

15 


CALL LNORMCSEED, Yl, 1, 1, 0) 



1 

16 


Yl = Yl * SIGMAY 



1 

17 


TOTTIM = 0. 0 



1 

18 


DETTIM = 10000.0 



1 

19 


CALL LRND(SEED, COURSE, 1, 1, 0) 



1 

20 


COURSE = COURSE * 6.283185 



1 

21 


CALL LRND(SEED, U, 1, 1, 0) 



1 

22 


IF (U .LE. 1.0/3. 0) THEN 


1 

1 

23 


SPEED = AVGSP*(24 + SQRT(3024*U + 16))/56 


1 

1 

24 


ELSE 


1 

1 

25 


SPEED = AVGSP*(60 - SQRT(1528 - 1512*U))/28 

1 

1 

26 


END IF 



1 

27 


CALL LEXPN(SEED, T, 1, 1, 0) 



1 

28 


LEGTM = T/LMBDA 



1 

29 


CALL LEXPN(SEED, AFTIM, 1, 1, 0) 



1 

30 

10 

CALL LNORMCSEED, ACFL, 1, 1, 0) 



1 

31 


ACFL = ACFL * 8.0 



1 

32 


FOM = FOM + ACFL 



1 

33 


IF (FOM .LE. 66.0 .OR. FOM .GE. 100.66) GO ' 

TO 10 


1 

34 

20 

CONTINUE 



1 

35 


READ(03,*) R, TL 



1 

36 


IF (FOM .GE. TL) THEN 


1 

1 

37 


Rl = R 


1 

1 

38 


TLl = TL 


1 

1 

39 


ELSE 


1 

1 

40 


DETRNG = Rl + (FOM - TL1)/(TL - TLl) 


1 

1 

41 


GO TO 30 


1 

1 

42 


END IF 



1 

43 


GO TO 20 



1 

44 

30 

CONTI.NUE 



1 

45 


REWIND 03 



1 

46 

35 

CONTINUE 



1 

47 


IF (AFTIM .LT. LEGTM) THEN 


1 

1 

48 


X2 = XI + SPEED * AFTIM * SIN(COURSE) 


1 

1 

49 


Y2 = Yl + SPEED * AFTIM * COS(COURSE) 


1 

1 

50 

40 

CONTINUE 


1 

1 

51 


READ(02,*) XC, YC 


1 

1 

52 


IF (XC .GT. 999.0) GO TO 50 


1 

1 

53 


CALL DTCT(DETRNG, XI, Yl, X2, Y2, XC, 

YC, DET, DSTDET) 

1 

1 

54 


IF (DET .EQ. 1) THEN 


2 

1 

55 


DT = TOTTIM + DSTDET/SPEED 


2 

1 

56 


IF (DT .LT. DETTIM) THEN 



39 




3 

1 

57 


DETTIM = DT 


3 

1 

58 


END IF 


2 

1 

59 


END IF 


1 

1 

60 


GO TO 40 


1 

1 

61 

50 

CONTINUE 


1 

1 

62 


REWIND 02 


1 

1 

63 


XI = X2 


1 

1 

64 


Yl = Y2 


1 

1 

65 


LEGTM = LEGTM - AFTIM 


1 

1 

66 


TOTTIM = TOTTIM + AFTIM 

} 

1 

1 

67 


CALL LEXPN(SEED, AFTIM, 1, 1, 0) 


1 

1 

68 

60 

CALL LNORMCSEED, ACFL, 1, 1, 0) 


1 

1 

69 


ACFL = ACFL * 8. 0 


1 

1 

70 


FOM = FOM + ACFL 


1 

1 

71 


IF (FOM .LE. 66.0 .OR. FOM .GE. 100.66) GO TO 60 

1 

1 

1 

72 

70 

CONTINUE 


1 

1 

73 


READ(03,*) R, TL 


1 

1 

74 


IF (FOM .GE. TL) THEN 


2 

1 

75 


R1 = R 


2 

1 

76 


TLl = TL 


2 

1 

77 


ELSE 


2 

1 

78 


DETRNG = R1 + (FOM - TL1)/(TL - TLl) 


2 

1 

79 


CO TO 80 


2 

1 

80 


END IF 


1 

1 

81 


GO TO 70 


1 

1 

82 

80 

CONTINUE 


1 

1 

83 


REWIND 03 


1 

1 

84 


ELSE 


1 

1 

85 


X2 = XI + SPEED * LEGTM * SIN(COURSE) 


1 

1 

86 


Y2 = Yl + SPEED * LEGTM * COS(COURSE) 


1 

1 

87 

90 

CONTINUE 


1 

1 

88 


READ(02,*) XC, YC 


1 

1 

89 


IF (XC . GT. 999.0) GO TO 100 

dsttSet; 

1 

1 

90 


CALL DTCT(DETRNG, XI, Yl, X2, Y2, XC, YC, DET, 

1 

1 

91 


IF (DET .EQ. 1) THEN 


2 

1 

92 


DT = TOTTIM + DSTDET/SPEED 


2 

1 

93 


IF (DT .LT. DETTIM) THEN 


3 

1 

94 


DETTIM = DT 


3 

1 

95 


END IF 


2 

1 

96 


END IF 


1 

1 

97 


GO TO 90 


1 

1 

98 

100 

CONTINUE 


1 

1 

99 


REWIND 02 


1 

1 

100 


XI = X2 


1 

1 

101 


Yl = Y2 


1 

1 

102 


AFTIM = AFTIM - LEGTM 


1 

1 

103 


TOTTIM = TOTTIM + LEGTM 


1 

1 

104 


CALL LRND(SEED, COURSE, 1, 1, 0) 


1 

1 

105 


COURSE = COURSE * 6.283185 


1 

1 

106 


CALL LRND(SEED, U, 1, 1, 0) 


1 

1 

107 


IF (U .LE. 1. 0/3.0) THEN 


2 

1 

108 


SPEED = AVGSP*(24 + SQRT(3024*U + 16))/56 


2 

1 

109 


ELSE 


2 

1 

110 


SPEED = AVGSP*(60 - SQRT(1528 - 1512*U))/28 


2 

1 

111 


END IF 

/ 

1 

1 

112 


CALL LEXPN(SEED, T, 1, 1, 0) 



40 


1 1 113 

1 1 114 

1 115 

1 116 
1 1 117 

1 1 118 

1 1 119 

1 120 
121 
122 


LEGTM = T/LMBDA 
END IF 

IF (TOTTIM .LT. SCHTM .AND. DETTIM . GT. 9999. 0) GO TO 35 
IF (DETTIM .LT. SCHTM) THEN 
NUMDET = NUMDET + 1 
WRITEdl,*) DETTIM 
END IF 
110 CONTINUE 

PRINDV, 'NUMBER OF DETECTIONS =', NUMDET 
STOP 


123 END 
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17 

18 

19 

20 
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22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 


SUBROUTINE DTCT(R, XI, Yl, X2, Y2, XC, YC, DET, DSTDET) 
REAL R, RSQ, XI, Yl, X2, Y2, XC, YC, XP, YP, M, DSQ 
REAL D12SQ, DIPSQ, D2PSQ, Z, DSTDET 
INTEGER DET 
RSQ = R * R 
IF (XI .EQ. 

XP = XI 
YP = YC 
ELSE IF (Yl 
XP = XC 
YP = Yl 
ELSE 

M = (Y2 


X2) THEN 


.EQ. Y2) THEN 


- XI) 
M*X1 ■ 


Y1)/(M + 1/M) 


X2) 

XP) 

X2) 


(Yl - Y2)*(Y1 - 
(Yl - YP)*(Y1 - 
(YP - Y2)*(YP - 
- Y1)*(YC - Yl) 


Y1)/(X2 
XP = (XC/M + YC + 

YP = M*XP - M*X1 + Yl 
END IF 

DSQ = (XC - XP)*(XC-XP) + (YC - YP)*(YC - YP) 
DET = 0 

IF (DSQ .LT. RSQ) THEN 
D12SQ = (XI - X2)*(X1 
DIPSQ = (XI - XP)*(X1 
D2PSQ = (XP " X2)*(XP 
IF ((XC - X1)*(XC - XI) + (YC 
DET = 1 
DSTDET = 0. 0 

ELSE IF ((DIPSQ . LT. D12SQ) 

DET = 1 

Z = SQRT(RSQ - DSQ) 

DSTDET = SQRT(DIPSQ) 

ELSE IF ((X2 - XC)*(X2 
^ THEN 

DET = 1 

Z = SQRT(RSQ - DSQ) 

DSTDET = SQRT(DIPSQ) 

END IF 
END IF 
RETURN 
END 


Y2) 

YP) 

Y2) 

.LT. 


RSQ) THE 


- Z 

XC) 


- Z 


AND. (D2PSQ .LT. D12SQ)) THEN 


+ (Y2 - YC)*(Y2 - YC) .LT. I^S(3 
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